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Abstract 


Signal  detection  in  a  large  segment  of  non- 
Gaussian  and  non-stationary  Arctic  under-ice  noise, 
which  contains  both  high  power  narrow-band  and 
impulsive  components,  is  examined.  It  is  shown  that 
the  correlation  functions  of  sub-segments  of  data 
change  significantly,  and  if  ignored,  can  degrade  the 
performance  of  a  detector.  For  a  false  alarm  proba¬ 
bility  of  0.05  and  a  known  constant  signal,  the 
matched  filter  was  on  average  17.6f?  better  than  a 
detector  designed  assuming  independent  noise  sam¬ 
ples.  It  is  also  shown  that  pre-processing  the  data 
with  an  adaptive  notch  filter,  then  using  the  matched 
filter  will  result  in  a  further  improvement  of  about 
6cc.  Additionally,  the  effect  two  different  signal 
shapes  have  on  the  performance  of  the  matched  filter 
is  examined.  )<>•/  ’  r  ' '  ‘  i  St 
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With  the  goal  of  creating  a  database  of  under-ice 
noise,  in  1980  a  multi-institutional  research  program  colla¬ 
borated  to  record  Arctic  under-ice  data.  The  data.  FR.A\! 
II.  was  recorded  on  23-2-1  April  1980.  from  a  pack  ice  floe, 
at  8G°.V  25°  IF.  An  omnidirectional  hydrophone,  radio 
linked  to  a  receiver,  w  as  suspended  to  a  depth  of  91m  in 
4000m  deep  water.  The  data,  recorded  on  an  analog  device, 
w  as  then  bandpass  filtered  from  0.01-5kHz.  The  particular 
data  segment  being  analyzed  w  as  recorded  on  the  23  April 
from  11:30  to  11:40  pnt.  Subsequent  to  the  bandpass  filter 
previously  mentioned,  the  data  was  passed  through  a  low 
pass  filter  with  cutoff  of  2.5  kHz.  and  digitized  with  a  sam¬ 
pling  frequency  of  10  kHz.  A  total  of  6.150.144  data  sam¬ 
ples  were  taken,  which  conveniently  breaks  down  to  6006 
records  with  each  record  containing  exactlv  1024  samples 
11.21- 

Several  authors  have  looked  at  this  particular  data 
segment  in  varying  degrees  of  detail.  Dwyer  [1]  took  several 
different  approaches  to  analyzing  the  under-ice  data.  He 
examined  the  first  four  central  moments  in  the  time 
domain  and  conclud"d  that  the  noise  is  non-stationary  and 
non-Gaussian.  A  similar  analysis  carried  out  in  the  fre¬ 
quency  domain  yielded  comparable  results.  The  energy 
spectrum  indicated  that  the  ice  noise  was  heavier  tailed 
than  a  Gaussian  source  processed  in  the  same  fashion,  in 
agreement  with  the  results  of  Greene  and  Buck  (3|.  Front 
an  examination  of  the  spectrum  and  the  time  domain  plot' 
of  the  data,  both  narrow-band  and  impulsive  component' 


were  identified.  The  existent  of  these  components  was 
not  surprising,  since  the  floe  was  moving  slowly  throughout 
the  experiment.  As  suggested  by  the  work  of  Milne  and 
Ganton  [4],  the  movement  of  the  floe  would  cause  rifting 
and  cracking  of  the  ice. 

Veitch  and  Wilks  [2]  took  a  slightly  different  approach 
since  their  objective  was  to  develop  a  statistical  model  for 
the  under-ice  noise.  Their  work  confirms  the  presence  of 
narrow- band  and  impulsive  components,  but  suggests  that 
the  background  noise  is  Gaussian.  By  examining  the  first 
four  central  moments  of  the  data  in  detail,  they  noted  a 
strong  correlation  between  records  with  large  kurtosis  and 
large  skew.  They  also  determined  that  a  large  kurtosis  was 
associated  with  high  amplitude  bursts.  In  addition,  it  was 
proposed  that  the  Arctic  data  could  be  modeled  in  the  fre¬ 
quency  domain  as  a  mixture  of  a  Gaussian  background  pro¬ 
cess.  a  sum  of  a  random  number  of  narrow-band  com¬ 
ponents.  each  with  random  phase  and  amplitude,  and  what 
is,  in  some  sense,  an  error  term.  The  main  drawbacks  of 
the  model  are  that  the  impulsive  component,  which  has 
been  verified  as  being  present,  is  not  modeled  specifically. 
Also  the  weights  and  number  of  sinusoidal'  components  are 
not  modeled  constructively:  that  is.  for  some  arbitrary  seg¬ 
ment  of  the  data,  the  number  of  tones  is  unknow  n. 

Recently.  Nielsen  and  Thomas  [5]  proposed  that  the 
univariate  statistics  of  the  Arctic  data  could  be  modeled  by 
the  generalized  Gaussian  family,  with  the  parameters  of  the 
model  estimated  from  the  data.  The  model  was  shown  to 
be  a  reasonably  accurate  representation  of  the  empirical 
probability  mass  functions  for  blocks  (5  records)  of  data, 
but  was  not  accurate  enough  for  the  resulting  Neyman- 
Pearson  optimal  detector  to  show  an  improvement  over  a 
simple  linear  structure. 

Time  porrelation 

Let  a  sample  function  of  a  continuous  parameter  ran¬ 
dom  process  be  sampled  at  the  Nvquist  frequency,  yielding 
samples  x,.  Define  the  time  correlation  functions  as 

1  v 
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Equivalently,  for  any  realization  x,  the  Fourier  transform  of 
the  realization  can  be  used  to  compute  the  power  spectrum, 
which  can  then  be  inverse  Fourier  transformed  to  obtain  an 
estimate  of  the  time  correlation  function.  The  Arctic  data 
was  sampled  at  10  kHz.  which  is  twice  the  Nvquist  fre¬ 
quency.  For  these  computation'  the  data  was  decimated  by 
a  factor  of  two  so  that  the  sampling  frequency  was  equal  to 
the  Nvquist  frequency:  thus,  the  correlation  functions  were 
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computed  using  .V=512. 

Since  time  correlation  functions  for  real  processes  are 
even  about  the  origin.  hm O,  the  graphs  of  typical  correla¬ 
tion  functions,  shown  in  Figures  1  -  4,  are  plotted  only  for 
<•» O.  1.  2 .  256.  In  general,  it  was  found  that  the  correla¬ 

tion  function  took  on  one  of  four  different  shapes.  The 
time  correlation  function  of  record  0001.  Figure  1,  shows  a 
mixture  of  a  decaying  low  frequency  component  with  a 
sir  'I  amplitude  high  frequency  component.  The  main 
variations  for  this  type  of  function  are  the  rate  of  decay 
and  the  amplitude  of  the  high  frequency  component. 


Figure  1  -  Time  Correlation  Record  0901 


Figure  2  shows  the  correlation  function  for  record 
090-5.  approximately  one-half  second  later  than  record  0901. 
The  correlation  drops  rapidly  as  k  increases,  and  then  has 
small  amplitude  variations  about  zero.  A  close  examination 
of  this  correlation  function  shows  that  significant  correla¬ 
tions  exist  for  the  first  three  shifts.  The  drastic  change  in 


Figure  2  -  Time  Correlation  Record  0905 


the  shape  of  the  correlation  function  in  a  short  period  of 
time,  as  illustrated  by  Figures  1  and  2.  is  not  unusual  for 
this  data. 

The  third  typical  correlation  function  is  shown  in  Fig¬ 
ure  3.  For  record  1205  the  correlation  is  a  mixture  of  a  low 
frequency  and  a  high  frequency  component.  The  result  is 
that  the  correlation  function  looks  like  a  corrupted  AM  sig¬ 
nal.  The  principal  variations  with  the  general  shape  are  in 
the  magnitude  of  the  correlation  for  small  k ,  and  the  decay 
rate  of  the  high  frequency  component. 
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Figure  3  -  Time  Correlation  Record  1205 

One  of  the  most  interesting  correlation  functions 
belongs  to  the  final  category.  In  Figure  4  the  correlation 
function  for  data  record  3905  is  shown.  The  correlation  for 
this  data  is  very  high,  and  periodic.  with  peak  correlations 
greater  than  0.5  for  shifts  greater  than  zero.  The  peak 
correlation  for  shifts  other  than  zero  were  found  to  be  plus 
or  minus  one  for  some  data  records,  and  the  envelope  shape 
varied  greatly.  The  main  factor  correlation  functions  of 
this  type  have  in  common  is  the  presence  of  a  single,  high 
powered  periodic  component. 

Adaptive  Notch  Filter 

The  correlation  between  data  samples,  which  has  been 
shown  to  be  significant  in  the  Arctic  under-ice  noise,  pro¬ 
vides  a  complicating  factor  in  the  problem  of  signal  detec¬ 
tion.  This  results  primarily  from  inability  to  model  the 
multivariate  density  if  the  data  is  known  to  be  non- 
Gaussian.  A  common  approach  to  working  with  highly 
correlated  noise  is  to  pre-whiten  the  data,  ideally  resulting 
in  independent  noise  samples.  Sub-optimal  pre-whiteners. 
which  are  designed  to  reduce  instead  of  eliminate  the  corre¬ 
lation  between  data  samples,  are  used  with  real  noise 
sources. 

For  environments  such  as  the  Arctic  data,  which  con¬ 
tain  a  combination  of  narrow-band  and  wideband  processes, 
an  adaptive  line  enhancer  [6j  is  commonly  used  to  separate 
the  two  processes.  Suppose  the  received  vector  Y  is  the 
sum  of  two  components;  a  white  background  component  N 
and  a  sinusoidal  component  S.  We  can  represent  this 
mathematically  as 
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Figure  4  -  Time  Correlation  Record  3005 
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In  Figure  5  we  design  the  filter  to  be  an  adaptive  notch 
filter  (ANF)  with  a  gain  of  unity  everywhere  except  at  a 
specified  band  of  frequencies,  where  the  gain  is  zero.  The 
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Figure  5  -  Adaptive  Line  Euhancer 


filter  output  will  be  the  white  noise  X  and  the  output  of 
the  difference  operator  will  be  the  sinusoid  5.  We  define 
the  ANF  as  follows  [7] 


N(z) 


z2  •+  bz  4-  1 
z 2  +  1  bz  +  ~f 


Then,  we  have 


bandwidth  2(1  — i),  -y « [0, 1  ] 

normalized  center  frequency  ==  cos- ‘(—6/2),  6  e[— 2,2j 

For  an  input  sequence  y*  and  output  sequence  of  the  ANF 
u>,  we  have 

wt  “  Vk  +  byk_ ,  +  y*_2  — 

For  a  fixed  notch  width,  the  only  adaptive  parameter  of  the 
filter  is  b.  which  controls  the  center  frequency  of  the  notch. 
To  derive  the  update  equation  for  6,  we  minimize  u>2  using 
a  gradient  search  algorithm,  then 

6t+ 1  “  <>k  +  step  tt>  (yt_,  -  -y  it?*,,) 

The  variable  step,  which  can  be  either  fixed  or  adaptively 
set.  controls  the  adaptation  speed  of  the  filter  and  the  vari¬ 
ance  of  the  filter  center  frequency  once  steady  state  has 
been  achieved.  That  is,  a  value  of  step  too  small  will  result 
in  a  long  time  for  the  filter  to  adjust  to  the  proper  center 
frequency,  however,  once  there,  only  small  changes  in  the 


filter  center  frequency  will  occur.  The  converse  results  will 
hold  for  a  value  of  step  too  large.  One  advantage  of  the 
ANF  is  that  it  can  be  implemented  easily  as  a  time  domain 
processer,  in  contrast  with  the  frequency  domain  pre¬ 
processor  proposed  by  Dwyer  [8j. 

As  an  example  of  the  effect  the  ANF  can  have  on  the 
correlation  function.  Figures  6-7  show  the  correlation 
function  for  record  -1201  before  and  after  filtering.  The 
ANF  essentially  removes  the  high  frequency  component 
which  dominates  the  correlation  function,  resulting  in  data 
which  is  much  less  correlated.  The  results  are  not  always 
as  dramatic  as  with  record  4201.  Specifically,  if  the  data 
contains  no  dominant  narrow-band  component,  little  differ¬ 
ence  in  the  correlation  function  after  filtering  is  observed. 
Also,  when  the  data  contains  multiple  narrow-band  com¬ 
ponents,  the  ANF  only  removes  the  highest  powered;  thus 
the  filtered  data  still  has  a  correlation  function  with 
periodic  components. 


Figure  6  -  Time  Correlation  Record  4201 


Figure  7  -  Time  Correlation  Filtered  Record  4201 


Signal  Detection 

Consider  the  following  simple  detection  problem.  Let 
A/  be  the  number  of  samples  taken  per  decision,  and  let  A' 
be  the  M-vector  of  samples  from  the  zero  mean  noise  pro¬ 
cess  and  s  be  the  known  signal  M-vector.  Let  j\f  represent 
the  M-variate  density  of  the  noise.  Define  a  binary 
hypothesis  problem  in  the  following  manner 


H0 :  A’ -A 

Hr.  A-A  +  s 

where  X  15  the  vector  received  by  the  detector  and  is  of 
length  A/.  Let 

Probability  (  decide  // ,  |  H o  true) 

3 «=  Probability  (  dectde  Hx  |  H |  true) 


The  Nevman-Pearson  criterion  can  be  written  as  follows. 
Fix  some  Qo-  Define  the  optimal  detector  to  be  the  detec¬ 
tor  which  maximizes  the  power  3  for  any  a  <  Oq.  Define 
the  likelihood  ratio  /  as  follows 


1(1) 


log 


fuiX-l) 

MX) 


Then  the  optimal  detector  can  be  written  as 

Hi)  >t-~hx 
Hi)  <  T  —*  H0 

The  threshold  T  is  chosen  to  satisfy  the  constraints  on  the 
power  3  and  on  the  level  a  [9]. 

.As  mentioned  before,  for  non-Gaussian  data  which  is 
correlated,  estimating  the  multivariate  density  of  the  noise 
in  ay  not  be  possible.  One  approach  taken  to  circumvent 
this  problem  is  to  assume  that  the  multivariate  statistics 
are  Gaussian:  then  using  the  Neyman-Pearson  criteria,  the 
optimal  likelihood  ratio  is  the  matched  filter.  If  R  is  the 
A/  x  A/  covariance  matrix  of  the  noise,  then  the  matched 
filter  is  defined  as 

Hi)  ~  str~' i 


For  independent  noise  samples,  £7V,-jV,  «>  6„<r  and  the 
covariance  matrix  R  —  <?  /,  where  /  is  an  M  x  M  identity 
matrix.  Consequently,  the  likelihood  ratio  reduces  to 
Hi)  *  The  resulting  detector  is  then  a  correlator. 

In  addition  only  the  variance  of  the  noise  <?  need  be 
estimated. 

For  correlated  data  the  estimation  of  the  elements  of 
the  covariance  matrix  is  necessary.  While  the  Arctic  noise 
does  not  appear  to  be  stationary  in  any  sense,  we  will 
assume  that  within  a  record  (that  is  1024  contiguous  sam¬ 
ples)  the  data  is  at  least  wide  sense  stationary'  and  ergodic. 
and  that  the  mean  and  variance  of  the  noise  are  known 
exactly.  LTnder  the  assumption  of  ergodicity,  the  correla¬ 
tion  function  is  equivalent  to  a  vector  of  expected  values. 
Pk  **  Define  the  covariance  matrix  of  the  noise 

as 


R, }  «  EX.Xj  -  Rj  ,  -  R,-j  —  Pi-j 

Each  record  contains  A’»»1024  samples;  thus  is  com¬ 
puted  for  1,  ...,  1023.  For  a  signal  length  of  A/,  only 
the  first  A/  correlations  are  used  in  the  formation  of  R. 
The  resulting  covariance  matrix  R  is  symmetric  and  Toe- 
plitz.  For  R  to  be  a  valid  invertible  cc,.  ..rianc"  matrix,  iw 


require  that  it  be  positive  definite. 

For  convenience,  the  data  was  normalized  such  that 
P0  =  1.  Thus,  the  covariance  matrix  is  o*R,  where  cr  is 
the  variance  of  the  record  and  R  is  the  normalized  covari¬ 
ance  matrix.  Henceforth,  we  refer  to  the  normalized  covari¬ 
ance  matrix  as  the  covariance  matrix.  A  signal  length  of 
M  —  4  was  used  for  all  simulations.  For  the  entire  data 
set.  6006  covariance  matrices  were  computed.  The 
minimum  eigenvalue  was  always  strictly  greater  than  zero, 
insuring  that  R  was  indeed  positive  definite,  and  typical 
values  for  the  four  eigenvalues  were  0.0077,  0.163,  1.216, 
and  2.581.  In  addition,  a  comparison  of  covariance 
matrices  showed  that  the  elements  changed  slowly  over 
time.  This  suggests  the  assumptions  of  wide  sense  stationar- 
ity  and  ergodicity  are  not  unreasonable  for  such  a  short 
segment  of  data  (1  record  »0.1  seconds). 

To  determine  the  effect  that  correlation  between  data 
samples  has  on  the  signal  detection  problem,  the  power  3 
for  the  following  detectors  was  compared  when  Af—4. 


[Dll 

D2 

D3 


Correlator  (R^cPl) 
Matched  Filter 
ANF  -  Matched  Filter 


For  [Dl]-[D3],  in  order  to  compensate  for  the  non¬ 
stationary  characteristic  of  the  noise,  parameters  were 
estimated  on  a  record-bv-record  basis.  For  [Dl],  o2  was 
estimated;  for  [D2]-[D3],  the  covariance  matrix  R  was 

estimated,  using  the  method  described  previously.  In  pari 
ticular,  for  ]D3],  the  ANF  was  used  as  a  preprocessor,  and 
the  covariance  matrix  was  estimated  using  the  filtered 
data.  For  convenience  the  data  was  grouped  into  blocks  of 
20  records.  Thus  each  block  contained  20480  data  samples, 
and  represented  a  time  interval  of  about  2  seconds.  We 
had  a  total  of  300  blocks,  and  the  performance  was 
evaluated  for  each  block. 


In  Table  1.  the  mean  and  variance  of  the  ratio  of  the 
power  levels  for  [Dl]  to  [D2j  are  shown,  along  with  the 
minimum  and  maximum  values  of  the  ratio,  and  the  perfor¬ 
mance  improvement.  We  defined  the  performance  improve¬ 
ment  to  be 


I{Dj,  />.)- 
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The  results  for  detectors  [D2]  and  [D3]  are  presented  in 
Table  2.  Examining  Tables  1  and  2.  we  se  that  the  largest 
improvements  in  performance  are  for  a  —  0.05;  and  while 
the  performance  improvement  decreases  as  a  increases,  the 
improvement  is  still  significant  for  a  »»  0.20.  For  a  —  0.05, 
[D2]  is  approximately  17.8%  better  and  [D3]  is  approxi¬ 
mately  29.2%  better,  on  average,  than  [Dl].  For  a  ■■  0.20. 
the  improvements  are  smaller,  but  at  12.4%  and  19.3%, 
they  are  still  significant. 


Given  that  the  Arctic  data  is  highly  correlated  at 
times,  the  improvement  that  the  matched  filter  has  over 
the  simple  correlator  is  understandable.  The  additional 
improvement  which  occurs  when  the  ANF  is  used  as  a  pre¬ 
processor  is  not  as  c^ar.  ?*  is  well  known  that  the  output 
of  a  linear  recursive  filter  is  more  Gaussian  than  the  input. 
Thus  using  the  ANF  as  a  pre-processor  not  only  whitens 
the  data,  by  removing  the  highest  powered  periodic  com- 
nonent  hut  al«o  retails  in  data  which  i«  more  Gaussian.  It 


Ratio  of  Power  Levels  between  [Dl]  and  [D2] 
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Table  1 


Ratio  of  Power  Levels  between  [D2j  and  [D3] 


Table  2 


is  clear  from  these  results,  that  disregarding  the  correlation 
between  data  samples  will  result  in  a  performance  degrada¬ 
tion. 

The  minimum  and  maximum  values  of  the  ratios  are 
also  interesting.  While  the  maximum  values  of  the  ratio 
are  near  unity,  indicating  that  the  detectors  have  similar 
performance,  the  minimum  values  range  from  0.58  to  0.75, 
indicating  that  under  some  conditions  one  detector  is  much 
better  than  the  other.  In  particular,  blocks  80-90  and 
170-210  can  be  shown  by  an  analysis  of  their  correlation 
functions  to  contain  periodic  components.  For  these  cases. 
[D2]  is  about  25rc  better  and  [D3]  is  approximately  47cy 
better  than  [Dl]  for  a  =  0.20.  Figures  8-9  are  plots  of 

/^D2\  and  02\  /®,DSi  foe  a  “0.20  to  illustrate  more 
clearly  the  improvements. 

One  may  challenge  the  validity  of  using  the  same  data 
to  estimate  the  covariance  matrix  and  compute  the  detec¬ 
tion  probabilities.  Tests  were  run  where  the  covariance 
matrix  was  updated  every  record,  every  fifth  record,  every 
tenth  record,  and  every  twentieth  record.  It  was  found 
that  the  performance  levels  did  not  change  appreciably  for 
updates  occurring  as  slowly  as  every  ten  records.  If  the 
update  was  performed  every  twenty  records,  once  a  block, 
only  a  small  performance  loss  occurred. 

For  all  computations,  a  constant  signal  was  used  such 
that  the  sum  of  the  elements  of  the  signal  vector  equaled 
the  noise  standard  deviation,  i.e.  input  S/N“l.  For  the 
matched  filter  [D2].  the  output  S/N  can  be  written  as 

5  T  1 

7?  -  *  *  * 

For  any  positive  definite  matrix  A  such  that  Xmm.  \nax 
are  the  minimum  and  maximum  eigenvalues  of  A,  b> 
Rayleigh's  principle  we  have  [10] 
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block  number  (20  records/block) 

Figure  8  -  Ratio  of  Power  Levels  between  pi)  and  p2] 
for  a  “  0.20 


block  number  (20  records/block) 

Fiftuie  0  -  Ratio  of  Power  Levels  between  (D2j  and  [D3] 
for  a  -  0.20 
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0.20 

mean 

0.200 

0.311 

0.399 

0.474 

variance 

1.4le-03 

1 ,92e-U3 

2.23e-03 

2.21  e-03 

minimum 

0.136 

0.244 

0.337 _ 

0.408 

maximum 

0.337 

0.469 

0.557 

0.625 

KD2.D41 _ 

399.0 _ 

221.75 - 

150.44 _ 

110-94 _ 

Table  3 


On  average,  the  minimum  and  maximum  eigenvalues  of  the 
covariance  matrix  were  0.00709  and  2.581.  Using 
Rayleigh's  principle,  the  average  S/N  will  be  between  0.387 
and  130.0-1.  For  a  constant  signal  chosen  such  that 
sTs_  (T  (recall  the  actual  covariance  matrix  is  a2  ft,  so 
the  variance  cancels),  the  average  output  S/N  ■*  0.610, 
with  an  overall  minimum  and  maximum  of  0.311  and 
2. -185.  While  there  is  considerable  variation  of  the  output 
S/N  values,  a  constant  signal  results  in  an  output  S/N 
which  is  near  the  minimum  possible  value.  Thus,  an 
improvement  in  the  performance  is  achieved  over  the  sim¬ 
ple  correlator  even  when  the  signal  choice  is  nearly  the 
worst  possible. 

It  is  well  known  that  the  performance  of  the  matched 
filter  is  very  sensitive  to  the  shape  of  signal  chosen  [ll].  To 
illustrate  more  clearly  the  effect  signal  choice  has  on  the 
performance  of  the  matched  filter,  the  performance  level 
for  the  matched  filter  was  compared  when  the  signal  was 
constant  (previous  case)  and  when  the  signal  had  a  sign 
change  in  one  element  for  U=M. 

[D2l  Matched  Filter  -  s.  /  "  (l.  1.  1.  1) 

VA/ 

[D4]  Matched  Filler  -  s  =  ^ j ^  (1.  1,  1,  —1) 

The  output  S/N  for  [D4]  was  11.68  on  average,  which  is 
much  higher  than  the  constant  signal  case.  Consequently, 
the  performance  level  for  [Dd)  is  much  higher  than  that  of 

[D2].  Table  3  is  similar  to  that  of  Tables  1  and  2,  however, 
it  is  for  the  same  detector  but  different  signals.  It  is  clear, 
that  the  signal  with  the  sign  change  is  a  better  signal  than 
the  constant  signal  with  regard  to  the  detection  problem 
posed. 

Conclusions 

We  have  examined  the  performance  of  four  different 
detectors  in  Arctic  under-ice  noise.  Of  the  first  three 
detectors,  all  which  use  a  constant  signal,  the  first,  which  is 
based  on  an  assumption  of  independent  noise  samples,  has 
been  shown  to  have  the  worst  performance.  In  addition, 
the  performance  of  the  matched  filter  can  be  improved  by 
either  pre-processing  the  data  with  an  ANF  (detector  three) 
or  by  changing  the  signal  shape  (detector  four).  In  any 
case,  for  the  Arctic  data,  it  is  clear  that  the  correlation 
between  data  samples  is  significant  and  should  be  used  in 
the  design  of  a  detector. 
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